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Abstract 

Based on a dataset acquired by the BABAR experiment running on and near the T(4S') resonance 
from 1999-2002, an upper hmit is set on the rate of D^-D^ mixing using the decay mode D*~^ 
D^n~^, followed by a semi-leptonic decay of the D^. Results are compared to previous BABAR 
analysis using hadronic decays. We also set limits on the flavor-changing neutral current decays 
]jO _j. g+g- (^+^^) and the lepton-flavor violation decays D'^ e^/x^. 
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1 Overview 



The BABAR experiment, which is documented in detail elsewhere[l], has since its start in 1999 not 
only given results on i?-physics but also a series of new results in charm physics. With a cc cross 
section of 1.3 nb at the T(AS) resonance compared to the cross section of around 1.1 nb for i? 
production there is in fact a higher prompt charm production than B production. 

Here we present two new results from BABAR. The first one is a search for mixing between the 
neutral D meson states in the semi-lcptonic decay channel [2] while the other is a search for rare 
lepton decays of the neutral D meson[3]. Both are processes that, if seen with the current statistics, 
would be clear signs of physics beyond the Standard Model (SM). 

2 D^-D^ mixing 

Charm mixing is characterised by the two parameters x = Am/T and y = AT/2T, where Am 
(AT) is the mass (width) difference between the the two neutral D mass eigenstates, and F is the 
average width. We define the overall time-integrated mixing rate as Rmix = {x"^ + y^)/2. 

Mixing between the neutral charm mesons is, within the SM, heavily suppressed by the GIM 
mechanism. The expected mixing rate through box and di-penguin diagrams is O{10~^ — lO"^'^) 
but enhancements involving non-pertubative effects are possible. For a recent review of predictions 
for both the SM rate and possible New Physics contributions see[4]. 

To search for mixing the production flavour of the D-meson is tagged from the charge of the 
pion in the decay D*~^ — > D^n^ and the decay flavour is tagged from the charge of the electron in 
the decay K~e^Ue- Charge conjugation is implied everywhere. The decay where the pion 

and the electron have opposite charge (called the wrong sign mode), can only proceed when the 

oscillates into a before its decay. The right sign mode where the pion and electron have the 
same charge is used as a normalisation mode. In an analysis where the efficiency for right sign and 
wrong sign decays are identical, Rmix is simply given as the time-integrated ratio of the two decay 
modes. 

The analysis is based on a sample of 87 fb^^ and uses the K~e~^Ue sample while ignoring 

the less pure muon sample. The mass difference AM between the partially reconstructed D*^ 
candidate and the partially reconstructed D° candidate is together with particle identification the 
main selection criterion to obtain a pure sample. 

Separate neural networks, with input parameters specifically describing the daughters and 
globally describing the rest of the event, are used to select signal events and reconstruct the 
momentum vector. The neural networks, combined with charged kaon and electron particle iden- 
tification, provide a relatively pure selection of unmixed signal events and give a resolution in AM 
of 2.2 MeV/c2. 

The time distribution of the right sign control sample follows a simple exponential convoluted 
with a resolution function R, while the wrong sign signal has the form 

rws{t) = e-'^t'^R, (1) 

where t is measured in units of the lifetime. 

An unbinned extended likelihood fit is performed on the 2-dimensional distribution of signal and 
background in the variables t and AM; first on the right sign sample and then, with the shared 
parameters between the two datasets constrained, on the wrong sign sample. In Figs. 1 and 2, 




Figure 1: AM (left) and decay time (right) pro- 
jections of fit (solid lines) to RS data (points): 
(top left) AM signal region — unmixed signal 
(above dashed line), background (dashed line); 
(bottom left) magnified vertical scale AM full 
fit region — unmixed signal (white), D'^ back- 
ground (light grey), background (dark grey), 
zero- lifetime background (black); (top right) de- 
cay time signal region — signal and background 
components as in bottom left plot; (bottom 
right) magnified vertical scale decay time full fit 
region — signal and background components as 
in bottom left plot of this figure. 




Figure 2: AM (left) and decay time (right) pro- 
jections of fit (solid lines) to WS data (points): 
(top left) AM signal region — mixed signal 
(above dashed line), background (dashed line); 
(bottom left) AM fuh fit region — back- 
ground (white), zero-lifetime background (dark 
grey), non-peaking background (intermedi- 
ate grey), peaking Z?+ background (light grey), 
mixed signal (black); (top right) decay time 
signal region — signal and background com- 
ponents as in bottom left plot; (bottom right) 
magnified vertical scale decay time full fit re- 
gion — signal and background components as 
in bottom left plot. 



projections of the fit can be seen overlaid on the data. The result is a wrong sign signal yield of 
114 lb 61 events. 

Systematic errors arise mainly from the assumptions related to the shape of signal and back- 
ground in the AM variable and when added in quadrature add up to 34% of the statistical error. 
Combining the wrong sign yield with the right sign yield we get the final result 

R^i^ = 0.0023 ± 0.0012(stat) ± 0.0004(syst) (2) 
R^i^ < 0.0042 at 90% CL. (3) 

Systematics are taken into account by scaling the log likelihood curve for the fit to the wrong sign 
yield with the systematic error added in quadrature (Vl + 0.34^ = 1.06). The upper limit was 
calculated assuming a flat prior for the number of wrong sign events to be positive. In Fig. 3, the 
result is compared to previous results. 
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Figure 3: The BABAR limit compared to other 
prehminary and pubhshed results using the 
semi-leptonic decay to search for mixing[5, 6]. 
The thin solid circle shows the sensitivity of the 
BABAR analysis if zero wrong sign signal events 
were seen. For the comparison to the BABAR re- 
sults from a mixing search with hadronic decays 
modes [7, 8] it is assumed that there is no CP vi- 
olation and that the strong phase difference is 
zero. 



Figure 4: The dilepton invariant mass distribu- 
tion for each of the decay modes. The dashed 
lines indicate the optimised signal mass win- 
dows. 



3 Flavour- changing neutral current and lepton-flavour violating 
decays 

In this analysis, a search is performed for the flavour-changing neutral current (FCNC) decays 
— > e~^e~ and /^^^~ and the lepton-flavour violating (LFV) decays — > e^fi^. In the 

SM, the FCNC decays are highly suppressed by the GIM mechanism and the LFV decays are strictly 
forbidden. Compared to rare decay searches in the K and B sector, rare D decays are sensitive 
to new physics involving the up-quark sector such as certain i2-parity violating super symmetric 
models [9]. 

As in the previous analysis the is required to originate from a D*~^, but this time to ensure 
as clean a sample as possible. For the same reason, the is required to have a momentum above 
2.4 GeV/c in the T(4S') centre-of-mass frame to reduce background from combinatorics involving 
the decay products of B mesons. Electrons (muons) are identified with an efficiency of 95% (60%) 
with a hadron misidentification probability of 0.2% (2%) as measured on a r decay control sample. 

The decay — > 7r"^7r~ is used as a control sample as it has very similar kinematics and, as 
such, the systematic errors can be minimised. Apart from the particle identification, the selection 
of the control channel, is identical to the criteria used for the signal. 

Based on a sample of 122 fb~^ and after optimisation of the selection criteria the events seen in 



Fig. 4 remain. The background is estimated from the sidebands with a looser selection applied and 
then scaling it to the final selection taking the small correlation between the criteria into account. 

We do not see any signal in any of the channels and the branching fraction upper limits have 
been calculated using an extension of the Fcldman-Cousins method [10] that avoids the unwanted 
effect of the Feldman- Cousins method[ll] that the UL for a search can go down in case of an 
upwards fluctuation in the expected background. Our result and a comparison to previous published 
results[12, 13] can be seen in Table 1. 



^ e+e- £>" ^ 1^+11- e^//^ 

002 3.34 ± 0.31 021 

N^^"^^ 2.21 ±038 1.28 ±032 1.93 ±0.36 

iVbg 2.23 ± 038 4.63 ±0.45 2.14 ±0.36 

iVobs 3 1 

UL at 90% CL 1.2 x 10'^ 1.3 x 10"^ 8.1 x 10"^ 

Previous published limit [12, 13] 6.2 x IQ-*^ 2.0 x 10"*^ 8.1 x 10~" 



Table 1: The summary of the number of expected background events {N^^g), number of observed 
events (A^obs)) and the branching fraction upper limits at the 90% confidence level for each decay 
modes. The uncertainties quoted here are total uncertainties. The uncertainty of N^^ is negligible 
for the ee and e/x decay modes. 
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